Solitons in systems of coupled scalar fields 
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We present a method to obtain soliton solutions to relativistic system of coupled scalar fields. This 
is done by examining the energy associated to static field configurations. In this case we derive a set 
of first-order differential equations that solve the equations of motion when the energy saturates its 
lower bound. To illustrate the general results, we investigate some systems described by polynomial 
interactions in the coupled fields. 
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Soliton solutions are known to be of great interest in 
several branches of physics. They usually apper as solu- 
tions of the equations of motion that describe a dynam- 
ical system containging nonlinear interactions among its 
degrees of freedom. The issue here is that the basic prin- 
ciples evoked to find the physical contents of the system, 
generally lead to a set of coupled nonlinear second-order 
differential equations of motion. Unfortunately, however, 
there are no general rules for investigating the presence 
of soliton solutions for such a set of differential equations. 

In this paper we present a method for investigating 
the presence of soliton solutions in relativistic systems of 
coupled scalar fields. The general system we are going 
to consider is important in quantum field theory, and in 
other branches of physics. In condensed matter, for in- 
stance, it may represent the continuum version of quasi- 
one-dimensional structures that require at least two de- 
grees os freedom to be described. 

To be explicit, let us consider a general field theory 
decribing a pair of interacting real scalar fields in bidi- 
mensional spacetime. This system is governed by the 
following Lagrangian density. 
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where U(<j>,x) is the potential, which identifies the par- 
ticular system one is interested in. We are using natural 
units and so the time t — x° = xq and the space coordi- 
nate x = x 1 = — x\ have the dimension of inverse energy, 
and the fields are dimensionless. 

The equations of motion that follow from Q are given 

by 
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For static field configurations <j) — 4>{x) and x — x( x ) we 
get 
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The potential U(cj>, x) m general is a nonlinear function of 
the coupled fields, and so the above equations of motion 
are coupled nonlinear second-order ordinary differential 
equations. In the standard way l3i iM °f search- 
ing for soliton solutions, one has to deal with the above 
equations of motion. 

Here, however, we shall follow another route In 
this case we investigate the energy corresponding to the 
static field configurations. For the above system, the 
energy can be written as 
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In order to present a general procedure, let us now in- 
troduce tow arbitrary but smooth functions F — F(<p, x) 
and G — G(<f>, x). Then, we work with the tow first terms 
in expression JSJ to rewrite the energy in the following 
form, 

E = E b + E', (7) 

where E\, is given by 

E b = |^(oo),x(oo)]-^(-oo), X (-oo)]| 

+ |g[^(oo),x(oo)]-g[0(-oo),x(-oo)]| (8) 

with x) = 4>F{4>, x) and 0(0, X ) = xG{4>, x), and E' 
has the form 
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From expression we recognize that if the potencial 
has the specific form 

TT 1/ A dF dG\ 2 \f„ A dF 9G x2 

u =2{ F + ^ + x oj) M G + *fl* +x a*, 

(10) 
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the energy is saturated to its bound Ef, given in (JSJ for 
static field configurations satisfying the following set of 
first-order equations 
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We now work with Eqs. (|llfl and (|12|) to verify that any 
pair (</>, x) that solves this set of first-order equations 
also solves the second-order equations of motion and 
(J3J) when the potencial is given by expression QlUp. This 
procedure then changes the task of solving sets of second- 
order differential equations to the task of solving sets of 
first-order differential equations. The price one pays is 
that this is not always possible, since the potencial must 
have the specific form given by expression (|ll)fl . When 
this is the case, however, we search for the minimum en- 
ergy configurations that will certainly play an important 
physical role, and whose energy one calculates immedi- 
ately via expression JSJ. 

The above precedure is completely general, and can be 
generalized to systems of three or more fields straightfor- 
wardly. However, to investigate the presence of topolog- 
ical solitons we impose the restrictions 
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Here we are guided by the fact that parity symmetry 
is necessary (but not sufficient) to give rise to topo- 
logical solitons. Now, from (|l(Jfl and l|13|) one sees 
that the functions F and G must have opposite par- 
ity in each of the fields: if F{(j),x) — ±F(— 4>,x)i t nen 
G(4>, x) = TG(-4>, x), and if F(4>, X ) = ±F(4>, - X ), then 
G(0,x) = TG(0,-x). 

To investigate particular systems, we must specify the 
functions F and G. We do this by using dimensionless 
(a,b,. . . ) real and positive and dimensional (A, /i, . . .) real 
parameters - dimensional parameters have the dimension 
of energy. First we note that if one chooses F — F(4>) and 
G = G(x), the two scalar fields decouple. As a simple 
illustration, let us consider F = F(<f) and G — 0. Here 
we get 



U{<t>)= 2 \ F 



^t + F + 4>— - 
dx dd> 



We now choose F in the form 
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to get 



U {<!>) = -\ 2 (<f-a 2 ), 
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and so the first-order equation becomes 

i 2 ) = 0. (18) 

This is the 4> 4 system. It has the well known solutions 

4>{x) = a tanh[Aa(x + x)). (19) 

Since we are considering a > 0, the sing of A specifies the 
particular solution: A > for the kink, and A < for the 
antikink. The corresponding energy is given by 
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Another simple illustration is obtained by choosing F 
F(4>) in the form 



F{<f>) 
In this case we get 

and 



^(^ 2 -a 2 ). 
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g + A^ 2 - a 2 ) = 0. 
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This is the <fi 6 system. It has [6j soliton solutions, which 
we can write as 

(j) 2 (x) = ia 2 {l + tanh[Aa 2 (x + x)]} . (24) 

Note that the sign of A identifies different solutions. The 
correponding energy is 



E=\\\\a\ 



(25) 



To go to a more general situation we now consider F — 
F((j>) and G = G(<fi,x)- From the reasoning just below 
Eqs. l(T3|l . we note that G(4>,x) — —G{(f>,— X )i and that 
F((j)) and G(</>, x) must have an opposite parity in 4>. 
Then, if one chooses for F((f) the same function we have 
already considered in Ijl6(l , the simplest choice for G(<f>, x) 
is given by 



G i<P,x) = ^M^X- 
Here the system is specified by potencial 



(26) 



A(0 2 - a 2 ) + ^ X 2 



^W. (27) 



The corresponding first-order equations that follow from 
(|llf> and (|12|) are given by 



- + X((j) -a) + -/* X 



0. 



dx 



+ = o. 
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The above system comprises soliton solutions. To find 
some of them, first one set = 0. This implies that 
X 2 = 2(A//i)a 2 , and so no soliton solution is found; note, 
however, that the parameters A and [i must have the same 
sign. Next, we set \ = 0. In this case, the first-order 
equations and reduce to the single equation 
(|i8|l . Here the pair of solutions is given by 

4>{x) — atanh[Aa(:r +3;)], x( x ) = 0- (30) 

The corresponding energy is calculated immediately, and 
give the same value (1201) we have already obtaind for the 
(f> system. 

To find other soliton solutions we use the trial orbit 
method introduced by Rajaraman Q to choose an orbit. 
Here, however, we are dealing with first-order equations 
and this seems to ease the task of choosing orbits. For 
instance, from Eq. I|29|l one sees that orbits like (j> 2 — a 2 — 



bx and 



bx + C X are good ones. Moreover, 



the orbit one chooses must be compatible with the first- 
order equation and J2HJ), and this provides a direct 
verification of the trial choice one is testing. We follow 
this reasoning to choose 

_ a 2) + = ^2 _ a 2 } (31) 



In this case Eq. (|28|l changes to 



dx 



+ AiW> -a ) =0 



(32) 



We solve this equation to get 

4>(x) — a tanh[^ta(a; + x)], (33) 

and so 

X(x) = ±a[2(A//i- l)] 1/2 sech[[ia(x + x)}, (34) 

which requires |A| > \fi\. To calculate the energy corre- 
sponding to the above solutions we use expression |(HJ): 
it is independent of fi and has the same value we have 
already calculated in (|2"U|> . 

To get another pair of solutions we set 



t> 2 - a 2 = x 2 ± 2 



ab 



X, 



(35) 



in which b > 1. In this case the first-order Eq. 1|29|) 
changes to 

l + "(° 2±2 W + ^) 1/2 ^ ' (36) 



which can be integrated to give 



X(x) = To 



b + cosh[(ia(x + x)] ' 



(37) 



and so we get 



(x) = a 



sinh[/iet(a 



b + cosh[(j,a(x + x)] 



(38) 



which requires fi = 2A. Once again, the energy is given 
by expression l|20|l. and now it is independent of b. 

As we have explicitly shown, the above soliton solu- 
tions J20I), and 1(23, and (123 and ^ have the 
same energy. This is so because they belong to the same 
topological sector. 

To investigate another system, let us choose F((f>) in 
the form give in (|21|) . In this case the simplest choice for 
G(<p, x) is given by 



G{<P,x) = t^0 2 X, 
and so the system is specified by the potential 



(39) 



U{4>, X ) = \[X^ 2 - a 2 )+^ X 2 ] 2 + ^W- (40) 

The corresponding first-order equations that now follow 
are given by 

£ + \^ 2 -a 2 )+^x 2 = 0, (41) 



— + ^x = o. 

dx 



(42) 



The above system comprises soliton solutions. To find 
some of them, first one sets <f> — 0. This implies that x 
is constant, and so no soliton solution is found. Next, we 
set x — 0, which reduces the above set of first-order equa- 
tions to the singles Eq. I|23f) . Here the pair of solutions 
is given by 



(j) 2 (x) = ^a 2 {l + tanh[Aa 2 (x + x)}}, 
X(x) = 0. 



(43) 



The correponding energy is calculated immediately, and 
gives the same value l|25|l we have already obtained for 
the 4> 6 system. 

To find other soliton solutions to (|41fl and (|42f> we fol- 
low the reasoning we have already considered in the for- 
mer system. Here we set 

A X (0 2 -a 2 )+^ X 2 =V<t>{4> 2 -a 2 ). (44) 
In this case the first-order equation (|41|l changes to 
dtp 



dx 



+ tMf>(<}> 2 - <f) = o. 



(45) 



The solutions to this equation are 

cj) 2 (x) = ia 2 {l + tanh[/ia 2 (x + x)]}, (46) 
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and so we get 



1 



r(z) = 2^(A///-l){l-tanh[/i^(x + ar)]}, (47) 

which requires A and /i to have the same sign, and |A| > 
To calculate the energy corresponding to the above 
solutions we use expression |J5}: it is independent of /i and 
has the same value we have already calculated in l|25|l . 

We now investigate one further system. Here we choose 
the function F = F((f>) yet as an odd function in <f>, but 
now we consider the simplest form 



(48) 



In this case the function G = G(<j>, x) has to be even in 
<j) and odd in x- As a simple choice we write 



G(0, X ) = ^A(0 2 -a 2 ) X . 



(49) 



In this case the system is specified by the following po- 
tential 

U{<j>, X ) = ^A 2 2 ( X 2 - a 2 ) 2 + \\\\<$? a 2 ) 2 , (50) 
and the corresponding first-order equations are given by 



^ + A(x 2 -a 2 )0 = O, 
^ + A(^ 2 -a 2 )x = 0. 



(51) 
(52) 



We now choose the orbit 2 = x 2 , and so the solutions 
are 

(f> 2 ( x ) = x 2 { x ) = - a 2 {l + tanh[Aa 2 (ir + x)]}. (53) 

Like in the former system, here the sign of A also iden- 
tifies different solutions. To calculate the corresponding 
energy, we use expression (JHJ to see that it has the same 
value we have already calculated in 1(25)1, 

An issue that now follows concerns investigating the 
stability of the static field configurations we have found. 
To give an explicit calculation, let us consider solutions 
(|3U[) to investigate the classical or linear stability in 
the simplest case. Here we write 



<j)(t, x) = 4>{x) + g(x) cos(eii), 
x(t, x) = x(x) + h(x) cos(e 2 t), 



(54) 
(55) 



where 4>(x) and x( x ) are given by l|3UI) . We now use the 
equation of motion J5J and to get 



cPh 
dz 2 



lb 2 



W2 A^ + f If + 1 ) S6Ch (Z) 



6sech 2 (z)] 5 = 0, (56) 



h = 0. (57) 



where z — \a(x + x) and uJ 2 — e 2 /X 2 a 2 , with i = 1,2. 
Eq. 1(561) is just the equation we get in the <f> system. 
It has |7j tow bound states, w 2 = and uj\ 1 — 3. Eq. 
(I57|l has bound states given by 



^l,n = ( 2 MA ~n)n, 



(58) 



where n = 0, 1, . . . < /i/A. There is no negative energy 
bound state, and so the set of solutions (|30|l is classically 
stable, at every point in the parameter space. This re- 
sult is interesting since one knows [8j that similar static 
field configurations that solve the system considered in 
Ref. 0] are linearly unstable, at least in some region in 
parameter space. 

For the set os solutions l(43|l . a similar investigation 
shows that it is also classically stable. The other solu- 
tions are harder to investigate, because the fluctuations 
become coupled and this greatly complicates the calcu- 
lations concerning classical stability. We shall return to 
this specifc issue in a separate work |9j. Here we just 
advance that all the soliton solutions we have found are 
classically stable. 

To end, we recall that we have investigated the en- 
ergy corresponding to static field configurations to offer 
a method to find soliton solutions to systems of coupled 
scalar fields. To illustrate the procedure, we have dealt 
with some specific systems, containing the forth and the 
sixth powers in the coupled fields. We then used the trial 
orbit method introduced by Rajaraman [j| to find soli- 
ton solutions. Evidently, the trial orbit method is not the 
only route. For instance, some very recent investigations 
jlfl lll| have introduced other methods, which can per- 
haps offer more soliton solutions for the system we have 
just investigated. 

As we have shown, the method seems to work well, 
at least for polynomial interactions between the cou- 
pled fields. However, further investigations are desired, 
mainly the one concerning the generalization to nonpoly- 
nomial interactions, which can be important to model 
periodic structures described by tow or more degrees of 
freedom. This and other related issues are now under 
consideration. 
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